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ABSTRACT 

We have measured resonance strengths and energies for dielectronic recombi- 
nation (DR) of Mg-like Fe XV forming Al-like Fe XIV via N = 3 N' = 3 core 
excitations in the electron-ion collision energy range 0-45 eV. All measurements 
were carried out using the heavy-ion Test Storage Ring at the Max Planck Insti- 
tute for Nuclear Physics in Heidelberg, Germany. We have also carried out new 
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multiconfiguration Breit-Pauli (MCBP) calculations using the AUTOSTRUC- 
TURE code. For electron-ion collision energies < 25 eV we find poor agreement 
between our experimental and theoretical resonance energies and strengths. From 
25 to 42 eV we find good agreement between the two for resonance energies. 
But in this energy range the theoretical resonance strengths are ~ 31% larger 
than the experimental results. This is larger than our estimated total experi- 
mental uncertainty in this energy range of ±26% (at a 90% confidence level). 
Above 42 eV the difference in the shape between the calculated and measured 
3s3p( 1 Pi)nZ DR series limit we attribute partly to the nl dependence of the de- 
tection probabilities of high Rydberg states in the experiment. We have used our 
measurements, supplemented by our AUTOSTRUCTURE calculations, to pro- 
duce a Maxwellian-averaged 3 — > 3 DR rate coefficient for Fe XV forming Fe XIV. 
The resulting rate coefficient is estimated to be accurate to better than ±29% (at 
a 90% confidence level) for k B T e > 1 eV. At temperatures of k-oT e ~ 2.5 — 15 eV, 
where Fe XV is predicted to form in photoionized plasmas, significant discrepan- 
cies are found between our experimentally-derived rate coefficient and previously 
published theoretical results. Our new MCBP plasma rate coefficient is 19 — 28% 
smaller than our experimental results over this temperature range. 

Subject headings: atomic data - atomic processes - plasmas - galaxies: active - 
galaxies: nuclei - X-rays: galaxies 



Introduction 



Recent Chandra and XMM Newton X-ray observations of active galactic nuclei (AGNs) 
have detected a new absorption feature in the 15-17 A wavelength range. This has been iden- 
tified as an unresolved transition array (UTA) due mainly to 2p — 3d inner shell absorption 
in iron ions with an open M- shell (Fe I - Fe XVI). UTAs h ave been obse rved in IRAS 



2002 



Kaspi et al 



178 flKaspi et al. 



2002; 



Behar et al. 



13349 +2438 flSako et alJ l200lh. Mrk-509 flPounds et al.ll200lh. NGC 3783 flBlustin et al. 



2003|), NGC 55 48 Jsteenbrugge et al.l[2003j), MR 2251^ 



20041) . I Zw 1 JGallo et all 120041 ). NGC 4051 (fPounds et all 120041 ) . and 



NGC 985 dKrongold et al.ll2005h . 



Based on atomic structure calculations and photoabsorbtion modeling, iBehar et al. 



(120011 ) have shown that the shape, central wavelength, and equivalent width of the UTA 
can be used to diagnose the properties of AGN warm absorbers. However, models which fit 
well absorption features from second and third row elements cannot reproduce correctly the 
observed UTAs due to the fourth row element iron. The models appear to predict too high 
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an ionization level for iron. iNetzer et al.l (120031 ) attributed this discrepancy to an underesti- 
mate of the low temperature dielect ronic re c ombin atio n (PR) rate coefficient s for Fe M-shell 
ions. To investigate this possibility iNetzerl (120041 ) and iKraemer et al.l (120041 ) arbitrarily in- 
creased the low temperature Fe M-shell DR rate coefficients. Th eir model res u lts ob tained 
with the modi fied DR rate coe fficients support the hypothesis of INetzer et al.l (120031 ). New 
calculations by iBadnelll (j2006al ) using a state-of-the -art theoretical method disscused in § [5] 
further support the hypotesis of INetzer et al.l (120031 ). 



Astrophysical mode l s cur rently use the DR data for Fe M-shell ions recommended 



Jacobs et al 



(11977f ) and 



by lArnaud fe Raymond! ( 19921) . Thes e data are based on theoretical DR calculations by 



Hahnl (I1989I ). The emphasis of this early theoretical work was 



on producing data for modeling collisional ionization equilibrium (sometimes also called 
coronal equilibrium). Under these conditions an ion forms at a temperature about an 
order of magnitude higher than the temper ature where it forms in ph otoionized plasmas 
( jKallman fe Bautistal 120011 ) . The use of the lArnaud fe Raymond! (119921 ) recommended DR 
data for modeling photoionized plasmas is thus questionable. Benchmarking by experiment 
is highly desirable. 

Reliable experiment ally- derived low temperature DR rate coefficients of M-shell iron 
ions are just now becoming available. Until recently, the only published Fe M-shell DR mea- 
surements were for Na-like Fe XVI (Linkemann et al. 1995; Miiller 1999; here and throughout 
we use the convention of identifying the recombination process by the initial charge state 
of the ion). The Na-like measur e ments were followed up wit h modern theoretical calcula- 



tions flGorczvca fe Badnell! Il996l : iGul 120041 : lAltun et al. 



2007). Additional M-shell experi- 



mental work also exis t s for Na-like Ni XVIII ( iFogle et al.ll2003l ) and Ar-like Sc IV and Ti V 
( iSchippers et al.lll998l . 120021 ) . We have undertaken to measure low temperature PR for other 
Fe M-shell ions . Our results for Al-like Fe XIV are presented in ISchmidt et al.l (120061 ) and 
Badnelll ( 12006bl ). The present paper is a continuation of this research. 



DR is a two-step recombination process that begins when a free electron approaches 
an ion, collisionally excites a bound electron of the ion and is simultaneously captured into 
a Rydberg level n. The electron excitation can be labeled NL —>■ N'l'-, where N is the 
principal quantum number of the core electron, I its orbital angular momentum, and j its 
total angular momentum. The intermediate state, formed by simultaneous excitation and 
capture, may autoionize. The DR process is complete when the intermediate state emits a 
photon which reduces the total energy of the recombined ion to below its ionization limit. 

In this paper we present experimental and theoretical results for AN=N' — N = DR 
of Mg-like Fe XV forming Al-like Fe XIV. In specific we have studied 3^3 DR via the 
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resonances: 



Fe i4+ (3 S 2 [ i S ]) + e 



' Fe 13+ (3s3p[ 3 P ° il2 ; 1 P 1 }nl) 
Fe 13+ (3s3d[ 3 Di, 2 , 3 ; 1 D 2 ]nl) 
Fe 13+ (3j9 2 [ 3 P ,i, 2 ; l D 2 - ^S Q ]nl) 

Fe 13+ (3p3 C /[ 3 D° 2i3 ; 3 F 2 ° 3i4 ; 3 P \ 2 ; l P°; ^F°]nl) 
^ Fe 13+ (3rf 2 [ 3 P ,i,2; 3 P 2 ,3,4; l D 2] l G A] l S ]nl) 



Possible contributions due to 3s3p 3 P metastable parent ions will be discussed below. Tabled] 
lists the excitation energies for the relevant Fe XV levels, relative to the ground state, that 
have been considered in our theoretical calculations. In our studies we have carried out 
measurements for electron-ion center-of-mass collision energies E cm between and 45 eV. 

Our work is motivated by the "formation zone" of Fe M-shell ions in photoionized gas. 
This zone may be defined as the temperatu re range where the frac tional abundance of a given 
ion is greater than 10% of its peak value (ISchippers et al.l I2004T ) . We adopt this definition 



for this paper. ISavin et al.l (119971 . Il999l . l2002al Jbl. 120061 ) defined this zone as the temperature 



range where the fractional abundance is greater t han 10% of the total elemental abundance. 
This is narrower than the ISchippers et al.l (120041) definition. Fo r Fe XV the wider definition 
corresponds to a k B T e rs 2.5-15 eV (IKallman fe Bautistall200ll ). It should be kept in mind 
that this temperature range depends on the accuracy of the underlying atomic data used to 
calculate the ionization balance. 

The paper is organized as follows: The experimental arrangement for our measure- 
ments is described in § [2j Possible contamination of our parent ion beam by metastable 
ions is discussed in § [31 Our laboratory results are presented in § HI In this section the 
experimentally-derived DR rate coefficient for a Maxwellian plasma is provided as well. 
Theoretical calculations which have been carried out for comparison with our experimental 
results are discussed in § El Comparison between the experimental and theoretical results is 
presented in § [61 A summary of our results is given in § [71 



2. Experimental Technique 



DR measurements were carried out at the heavy-ion test storage ring (TSR) of the 
Max-Planck Institute for Nuclear Physics (MPI-K) in Heidelberg, Germany. A merged 
beams technique was used. A beam of 56 Fe 14+ with an energy of 156 MeV was provided 
by the MPI-K accelerator facilit y Ions were inject ed into the ring and their energy spread 
reduced using electron cooling (IKilgus et al.lll990l ). Typical waiting times after injection 
and before measurement were ~ 1 s. Mean stored ion currents were ~ 10 fiA. Details of 
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the experimental setup have been given elsewhere (IKilgus et al.l Il992l ; lLampert et al.l 11996 
Schippers et~al~1ll99sl . I2OO0I boOlh . 



Recently a sec ond electron beam has been installed at the TSR (ISprenger et al.l 12004 



Kreckel et al.ll2005l ). This allows one to use the first electron beam for continuous cooling 
of the stored ions and to use the second electron beam as a target for the stored ions. In 
this way a low velocity and spatial spread of the ions can be maintained throughout the 
course of a DR measurement. The combination of an electron cooler and an electron target 
can be used to scan energy-dependent electron-ion collision cross sections with exceptional 
energy resolution. In comparison to the electron cooler, the electron source and the electron 
beam are considerably smaller and additional procedures, such as the stabilization of the 
beam positions during energy scans and electron beam profile measurements, are required 
to control the absolute luminosity product between the ion and electron beam on the same 
precise level as reached at the cooler. The target electron beam current was ~ 3 mA. The 
beam was adiabatically expanded from a diameter of 1.6 mm at the cathode to 7.5 mm 
in the interaction region using an expansion factor of 22. This was achieved by lowering 
the guiding magnetic field from 1.28 T at the cathode to 0.058 T in the interaction region 
thus reducing the transverse temperature to approximately 6 meV. The relative electron- 
ion collision energy can be precisely controlled and the recombination signal measured as a 
function of this energy. We estimate that the uncertainty of our scale for E cm is < 0.5%. 

The electrons are merged and demerged with the ion beam using toroidal magnets. 
After demerging, the primary and recombined ion beams pass through two correction dipole 
magnets and continue into a bending dipole magnet. Recombined ions are bent less strongly 
than the primary ion beam and they are directed onto a particle detector used in single 
particle counting mode. Some of the recombined ions can be field-ionized by motional 
electric fields between the electron target and the detector and thus are not detected. Here 
we assumed a sharp field ionization cutoff and estimated for Fe XV that only electrons 
captured into n max J$ 80 are detected by our experimental arrangement. 

The experimental energy distribution can be described as a flattened Maxwellian dis- 
tribution. It is characterized by the transversal and longitudinal temperatures Tj_ and Tji, 
respectively. The experimental energy spread depends on the electron-ion collision energy 
and can be approximat ed according to the formula AE = ([ln(2)A;BTj_] 2 -|-16 ln(2)i? cm A;B^||) 1 ^ 2 
(IPastuszka et al.lll996l ). For the comparison of our theoretical calculations with our experi- 
mental data we co nvolute the theoretic al results described in §[5] with the velocity distribution 
function given by iDittner et al.l (119861 ) to simulate the experimental energy spread. 



With the new combination of an electron target and an electron cooler we obtain in the 
present experiment electron temperatures of ksT± w 6 meV and ksTu « 0.05 meV. In order 
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to verify the absolute calibration of the absolute rate coefficient scale we also performed 
a measurement with the electron cooler using the previous standard method (Kilgus et al. 
1992, Lampert et al. 1996). We find consistent rate coefficients and spectral shapes, while 
the electron temperatures were larger by a factor of about 2 with the electron cooler alone. 
Moreover, because of the large density of resonances found in certain regions of the Fe XV 
DR spectrum the determination of the background level for the DR signal was considerably 
more reliable in the higher resolution electron target data than in the lower resolution cooler 
data. Hence, we performed the detailed analysis presented below on the electron target data 
only. 



D etails of th e experimental and data reduction procedure are given in ISchippers et al. 



( 120011 . 12004J ) and I Savin et al.l (120031 ) and reference therein. The baseline experimental un- 
certainty (systematic an d statistical) of the P R measurements is estimated to be ±25% at 
a 90% confidence level (ILampert et al.lll996l ). The major sources of uncertainties include 
the electron beam density determination, ion current measurements, and corrections for the 
merging and demerging of the two beams. Additional uncertainties discussed below result in 
a higher total experimental uncertainty as is explained in §§ 3 and 4. Unless stated otherwise 
all uncertainties in this paper are cited at an estimated 90% confidence level. 



3. Metastable Ions 



For Mg-like ions with zero nuclear spin (such as 56 Fe), the ls 2 2s 2 2p 6 3s3p 3 Pq level is 
forbidden to decay to the ground state via a one-photon transition and the multiphoton 
transitio n rate is negligible. Hence this level c an be considered as having a nearly infinite 
lifetime (IMarques. Parente. fe Indelicate! 1 19931 ; iBrage et al.l 119981 ) . It is possible that these 
metastables are present in the ion beam used for the present measurements. 

We estimate that the largest possible metastable 3 Pq fraction in our stored beam is 11%. 
This assumes that 100% of the initial Fe 14+ ions are in 3 Pj levels and that the levels are 
statistically populated. We expect that the J = 1 and 2 levels will radiatively decay to the 
ground state during the ~ 1 s between injection and measurement. The lifetimes of the 3 Pi 
and 3 P 2 levels are ~ 1.4 x 10" 10 s (Marques et al. 1993) and ~ 0.3 s (Brage et al. 1998), 
respectively. These decays leave l/9th or 11% of the stored ions in the 3 Po level. 

Our estimate is only slightly higher than the inferred metastable fraction for the ion 
beam used for DR measurements of the analogous Be-like Fe 22+ (Savin et al. 2006). The 
Be-like system has a metastable ls 2 2s2p 3 Pq state and following the above logic the stored 
Be-like ion beam had an estimated maximum 11% 3 Pq fraction. Fortunately, for the Be-like 



6 



measurements we were able to identify DR resonances due to the 3 Po parent ion and use the 
ratio of the experimental to theoretical resonance strengths to infer the 3 Po fraction. There 
we determined a metastable fraction of 7% ±2%. A similar fraction was inferred for DR 



measurements with Be-like Ti 18+ ions (ISchippers et al.l 120071 ). 



Using theory as a guide, we have searched our Mg-like data fruitlessly for clearly iden- 
tifiable DR resonances due to metastable 3 Po parent ions. First, followi ng our work in th e 
analogous Be-like Fe 22+ with its 2s2p 3 Po — * 2p 2 core excitation channel (jSavin et al.l 120061 ) . 
we searched for Fe 14+ resonances associated with the relevant 3s3p 3 Po — > 3p 2 core exci- 
tations. However, most of these yield only very small DR cross sections as they strongly 
autoionize into the 3s3p 3 Pj=i,2 continuum channels. These are energetically open at E cm 
greater than 0.713 eV and 2.468 eV, respectively (Tabled]). Hence, above E cm « 0.713 eV 
there are no predicted significant DR resonances for metastable Fe via 3s3p 3 P -> 3p 2 
core excitations. Below this energy the agreement between theory and experiment is ex- 
tremely poor (as can be seen in Fig. 1) and we are unable to assign unambiguously any DR 
resonance to either the ground state or metastable parent ion. Second, we searched for reso- 
nances associated with 3s3p 3 Po — > 3s3p 1 Pi, 3s3p 3 Po — ► 3s3p 3 Pi, and 3s3p 3 Po — > 3s3p 3 P2 
core excitation which are energetically possible for capture into the n > 14, 62, and 33 levels, 
respectively, and which may contribute to the observed resonance structures. The analogous 
2s2p 3 Pp — > 2s2y 1 P| and 2s2p 3 Pq — > 2s2p 3 P2 core excitations were seen for Be-like Ti 18+ 



(ISchippers et al.l 120071 ). However, again the complexity of the Fe XV DR resonance spec- 
trum (cf., Fig. 1) prevented unambiguous identification for DR via any of these three core 
excitations. Hence despite these two approaches, we have been unable to directly determine 
the metastable fraction of our Fe 14+ beam. 

Clearly our assumption that the 3 Pj levels are statistically populated is questionable. 
Ion beam generation using beam foil techniques are known to produce excited levels. The 
subsequent cascade relaxa t ion could potential ly populated the J levels non-statistically 
( Martinson fe Gaupp 1974 ; Quinet et ali 1999). Additionally the magnetic sub levels mj 
can be populated non-statistically ([Martinson fe Gaupplll974l ) which may affect the J lev- 
els. However, our argument in the above paragraphs that the 3 Pj levels are statistically 
populated yields 3 Po fractions of the analgous Be-like Ti 18+ and Fe 22+ of 11% while our 
measurements found metastable fractions of ~ 7% for those two beams. From this we con- 
clude either (a) that if 100% of the initial ions are in the 3 Pj levels, then the J = level is 
statistically under-populated or (b) that the fraction of initial ions in the 3 Pj levels is less 
than 100% by a quantity large enough that any non-statistical populating of the various J 
levels still yields only a 7% 3 Po metastable fraction of the ion beam. Thus we believe that 
our assumption provides a reasonable upper limit to the metastable fraction of the Fe + 
beam. 
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Based on our estimates above and the Be-like results we have assumed that 6%±6% of 
the Fe 14+ ions are in the 3s3p 3 Po metastable state and the remaining fraction in the 3s 2 1 S'o 
ground state. Here, we treat this possible 6% systematic error as a stochastic uncertainty 
and add it in quadrature with the 25% uncertainty discussed above. 



4. Experimental Results 

Our measured 3^3 DR resonance spectrum for Fe XV is shown in Figs. 1-8. The 
data (cry) represent the summed DR and radiative recombination (RR) cross sections times 
the relative velocity convolved with the energy spread of the experiment, i.e., a merged beam 
recombination rate coefficient (MBRRC). 

The strongest DR resonance series corresponds to 3s 2 1 Sq — > 3s3p l P\ core excitations. 
Other observed features in the DR resonance spectrum are possibly due to double core 
excitations discussed in §[TJ T rielectronic recombin at ion (TR), as this has been named, has 



been observed in Be-like ions (jSchnell et al.l l2003aU bl; iFogle et al.l 120051 ) . These ions are the 



second row analog to third row Mg-like ions. However in our data unambiguous assignment 
of possible candidates for the TR resonances could not be made. 

Extracted resonance energies Ei and resonance strengths Si for E cm < 0.95 eV are listed 
in Table [2] along with their fitting errors. These data were derived following the method 
outlined in Kilgus et al. (1992). Most of these resonances were not seen in any of the 
theoretical calculations for either ground state or metastable Fe 14+ . Hence their parentage 
is uncertain. The implications of this are discussed below. 

Difficulties in determining the non-resonant background level of the data contributed 
an uncertainty to the extracted DR resonance strengths. For the strongest peaks this was 
on the order of ~ 10% for E cm < 5 eV and ~ 3% for E cm > 5 eV. Taking into account the 
25% and 6% uncertainties discussed in §§ 2 and 3, respectively, this results in an estimated 
total experimental uncertainty for extracted DR resonance strengths of ±28% below ~ 5 eV 
and ±26% above. 

Due to the energy spread of the electron beam, resonances below E cm < ksT± cannot be 
resolved from the near eV RR signal. Here this limit corresponds to ~ 6 meV. But we can 
infer the absence of resonances lying below the lowest resolved resonance at 6.74 meV. For 
E cm ^$ &b2]|, a factor of up to ~ 2 — 3 enhanced MBRRC is observed in merged electron-ion 
beam experiments (see e.g., Gwinner et al. 2000; Heerlein et al. 2002). Here this temperature 
limit corresponds to E cm < 0.05 meV. As shown in Fig. [9[ at an energy 0.005 meV our 
MBRRC is a factor of 2.5 times larger than the fit to our data using the RR cross section 
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from semi-classical RR theory with quantum mechanical corrections (jSchippers et al.ll200ll ) 
and the extracted DR resonance strengths and energies. This enhancement is comparable 
to that found for systems with no unresolved DR resonances near eV (e.g., Savin et 
al. 2003 and Schippers at al. 2004). Hence, we infer that there are no additional significant 
unresolved DR resonances below 6.74 meV. Recent possible exp l anati ons for the cause of the 
enhancement near eV have been given by iHorndl et al.l (120051 . 120061 ) and reference therein. 



We have generated an experiment ally- derived rate coefficient for 3 —>■ 3 DR of Fe XV 
forming Fe XIV in a plasma with a Maxwellian electron energy distribution (Fig. [TO]) . For 
E cm < 0.95 eV we have used our extracted resonance strengths listed in Table [2j For 
energies E cm > 0.95 eV we have numerically integrated our MBRRC data after subtracting 
out the no n-resonant ba ckground. The rate coefficient was calculated us i ng th e methodology 
outlined in ISavinl ( 1999 ) for resonance strengths and in Schippers et al. ( 2001 ) for numerical 
integration. 

In the present experiment only DR involving capture into Rydberg levels with quantum 



numbers 



< 



80 contribute to the measured MBRRC. In order to generate a total AiV=0 



plasma rate coefficient we have used AUTOSTRUCTURE calculations (see § E]) to account 
for DR into higher n levels. As is discussed in more detail in § [BJ between 25-42 eV we find 
good agreement between the experimental and AUTOSTRUCTURE resonance energies. 
However, the theoretical results lie a factor of 1.31 above the measurement. To account 
for DR into n > n max = 80, above 42 eV we replaced the experimental data with the 
AUTOSTRUCTURE results (n max =1000) reduced by a factor of 1.31. Our resulting rate 
coefficient is shown in Fig. [TUJ 

Including the DR contribution due to capture into n > 80 increases our experimentally- 
derived DR plasma rate coefficient by < 1% for ksT e < 7 eV, by < 2.5% at 10 eV and by 
< 7% at 15 eV. This contribution increases to 20% at 40 eV, rises to 27% at 100 eV and 
saturates at « 35% at 1000 eV. Thus we see that accounting for DR into n > n max = 80 
levels has only a small effect a t temperatures of k^Tp. ~ 2 .5-15 eV where Fe XV is predicted 



to form in photoionized gas (IKallman fe Bautistal 120011 ). Also, any uncertainties in this 
theoretical addition, even if relatively large, would still have a rather small effect at these 
temperatures on our derived DR total rate coefficient. Hence, we have not included this in 
our determination below of the total experimental uncertainty for the experimentally-derived 
plasma rate coefficient at ksT e > 1 eV. 



The two lowest-energy resonances in the experimental spectrum occur at energies of 
6.74 meV and 9.80 meV with resonance strengths of 1.89 x 10~ 16 cm 2 eV and 1.01 x 
10~ 17 cm 2 eV, respectively (see Table 1 and Fig. [9]). As already mentioned, the parent- 
age for the two lowest energy resonances is uncertain. These resonances dominate the DR 
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rate coefficient for k B T e < 0.24 eV. The contribution is 50% at 0.24 eV, 16% at 0.5 eV, 
6.5% at 1 eV, 2.4% at 2.5 eV, and < 0.31% above 15 eV. At temperatures where Fe XV is 
predicted to form in photoionized plasmas, contributions due to these two resonances are 
insignificant. Because of this, we do not include the effects of these two resonances when 
calculating below the total experimental uncertainty for the experimentally-derived plasma 
rate coefficient at ksT e > 1 eV. 

An additional source of uncertainty in our results is due to possible contamination of 
the Fe XV beam by metastable 3 Pq ions. Because we cannot unambiguously identify DR 
resonances due to metastable parent ions, we cannot directly subtract out any contributions 
they may make to our experimentally-derived rate coefficient. Instead we have used our 
AUTOSTRUCTURE calculations for the metastable parent ion as a guide, multiplied them 
by 0.06 on the basis of the estimated (6 ± 6)% metastable content. We then integrated 
them to produce a Maxwellian rate coefficient and compared the results to our experimental 
results, leaving out the two lowest measured resonances at 6.74 and 9.80 meV. As discussed 
in the paragraph above, these two resonaces were left out because of the uncertainty in their 
parentage and their small to insignificant effects above 1 eV. The metastable theoretical 
results are 9.5% of this experimentally-derived rate coefficient at ksT e = 1 eV, 4.9% at 
2.5 eV, 2.2% at 5 eV, 1% at 10 eV and < 0.77% above 15 eV. 

In reality these are probably lower limits for the unsubtracted metastable contributions 
to our experimentally-derived rate coefficient. However, these limits appear to be reasonable 
estimates even taking into account the uncertainty in the exact value of the contributions 
due to metastable ions. For example, if we assume that we have the estimated maximum 
metastable fraction of 11%, then our experimentally-derived rate coefficients would have to 
reduced by only 9.0% at 2.5 eV, 4.0% at 5 eV, 1.8% at 10 eV, and less than 1.4% above 
15 eV. Alternatively, it is likely that theory underestimates the resonance strength for the 
metastable parent ions similar to the case for ground state parent ions (cf., Fig. 1). However, 
if the metastable fraction is 6% and the resonance contributions are a factor of 2 higher, then 
our experimentally-derived rate coefficients would have to reduced by only 9.8% at 2.5 eV, 
4.4% at 5 eV, 2.0% at 10 eV, and less than 1.5% above 15 eV. These are small and not very 
significant corrections. We consider it extremely unlikely that we have underestimated by a 
factor of nearly 2 both the metastable fraction and the metastable resonance contribution. 
Thus we expect contamination due to metastable 3 Po ions to have a small to insignificant 
effect on our derived rate coefficient at temperatures where Fe XV is predicted to form in 
photoinoized gas. 

Taking into account the baseline experimental uncertainty of 25%, the metastable frac- 
tion uncertainty of 6%, and the nonresonant background uncertainty of 10%/3%, all dis- 
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cussed above, as well as the uncertainty due to the possible unsubtracted metastable res- 
onances, the estimated uncertainty in the absolute magnitude of our total experimentally- 
derived Maxwellian rate coefficient ranges between 26% and 29% for ksTe > 1 eV. Here 
we conservatively take the total experimental uncertainty to be ±29%. This uncertainty 
increases rapidly below 1 eV due to the ambiguity of the parentage for the two lowest energy 
resonances and possible resonance contributions from metastable Fe XV which we have not 
been able to subtract out. 

We have fitted our experimentally-derived rate coefficient plus the theoretical estimate 
for capture into n > 80 using the simple fitting formula 

a DR (T e )=T-^J2 c ^ El/kBTe (2) 

i 

where q is the resonance strength for the ith fitting component and E{ the corresponding 
energy parameter. Table [3] lists the best-fit values for the fit parameters. All fits to the total 
experimentally-derived Maxwellian-averaged DR rate coefficient show deviations of less than 
1.5% for the temperature range 0.001 < ksT e < 10000 eV. 

In Table 3, the Experiment (I) column gives a detailed set of fitting parameters where 
the first 30 values of q and their corresponding Ei values are for all the resolved resonances 
for E cm < 0.95 eV given in Table 2. The parentage for these resonances are uncertain, though 
the majority are most likely due to ground state and not metastable Fe 14+ . It is our hope that 
future theoretical advances will allow one to determine which resonances are due to ground 
state ions and which are due to metastables. Listing the resonances as we have will allow 
future researchers to readily exclude those resonances which have been determined to be due 
to the metastable parent. The remaining 6 fitting parameters yield the rate coefficient due 
to all resonances for E cm between 0.95 and the 3s3p( 1 Pi)n/ series limit at 43.63 eV. In the 
Experiment (II) column of Table 3, the first six sets of and Ei give the fitting parameters 
for the first six resonances. The remaining sets of fit parameters are due to all resonances 
between 0.1 eV and the series limit. 



5. Theory 



The onl y published theoret ical DR rate coefficie nt for Fe XV which we are aware of is 



the work of I Jacobs et al.l (119771). Using th e work of iHahnl (119891 ). Arnaud and Raymond 



(1992) modified the results of iJacobs et al.l (119771 ) to take into account contributions from 
2p — 3d inner-shell transitions. The resulting rate coefficient of Arnaud and Raymond (1992) 
is widely used throughout the astrophysics community. 
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We have carried out new calculations using a state-of-the-art multiconfiguration Breit- 
Pauli (MCBP ) theor etical method. Details of the MCBP calculations have been reported in 
Badnell et al.l (120031 ). Briefly, the AUTOSTRUCTURE code was used to calculate energy 
levels as well as radiative and autoionization rates in the intermediate-coupling approxi- 
mation. These must be post-processed to obtain the final state level-resolved and total 
dielectronic recombination data. The resonances are calculat ed in the independent process 
and isolated resonance approximation ( ISeaton fa Storey! Il976l ). 



The ionic thresholds were shifted to known spectroscopic values for the 3 — > 3 transi- 
tions. Radiative transitions between autoionizing states were accounted for in the calculation. 
The DR cross section was approximated by the sum of Lorentzian profiles for all included 
resonances. The AUTOSTRUCTURE calculations were performed with explicit n values up 
to 80 in order to compare closely with experiment. The resulting MBRRC is presented for 
3 — > 3 core excitations in Figs. 1-8. 

The theoretical 3 — > 3 DR plasma rate coefficient was obtained by convolving calculated 
DR cross section times the relative electron-ion velocity with a Maxwellian electron energy 
distribution. Cross section calculations were carried out up to n max = 1000. The resulting 
Maxwellian plasma rate coefficient is given in Fig. [TUJ 

We have fit our theoretical 3^3 MCBP Maxwellian DR rate coefficients using Eq. [2j 
The resulting fit parameters are presented in Table [3j The accuracy of the MCBP fit is 
better than 0.5% for the temperature range 0.1 < ksTe < 10000 eV. This lower limit 
represents the range over which rate coefficient data were calculated. Data are not presented 
below (10 1 2: 2 )/11605 eV, w hich is estimated to be the lower limit of the reliability for the 
calculations (IBadnelll 120071 ). Here z = 14 and this limit is 0.17 eV. 



6. Discussion 



6.1. Resonance Structure 



As we have already noted, we find poor agreement between our experimental and the- 
oretical resonance energies and strengths for electron-ion collision energies below 25 eV. 
Theory does not correctly predict the strength of many DR resonances which are seen in 
the measurement. A similar extensive degree of disparity bet ween the theoretical and the 
measu red resonances was also seen in our recent Fe 13+ results ( jSchmidt et al.ll2006l : iBadnell 
2006bh . 



Some of the weaker peaks in our data below 1 eV may be due to the possible presence 
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of metastable Fe 14+ in our beam. But the estimated small metastable contamination seems 
unlikely to be able to account in this range for many of the strong resonances which are not 
seen in the present theory. Above ~ 1 eV, we expect no significant DR resonances due to 
metastable Fe 14+ (as is discussed in §[3]). 

In the energy range from 1 — 25 eV, the differences between experiment and theory are 
extensive. The reader can readily see from Figs. 1-8 that theory does not correctly predict 
the strength of many resonances which are observed in the experiment. This conclusion 
takes into account the by-eye shifting of the theoretical resonances energies to try to match 
up theory with the measured resonances. 

Between 25 — 42 eV we find good agreement between the experiment and theory for 
resonance energies. The AUTO STRUCTURE code reproduces well the more regular res- 
onance energy structure of high-n Rydberg resonances approaching the 3s3p( 1 Pi)n/ series 
limit. However the AUTOSTRUCTURE cross section lies ~ 31% above the measurements. 
This discrepancy is larger than the estimated ±26% total experiment al uncertainty in this 



energy range. A similar discrepancy with theory was found for Fe + (]Badnellll2006bh 



Theory and experiment diverge above 42 eV and approaching the 3s3p( 1 Pi)n/ se- 
ries limit. We attribute the difference in the shape between the calculated and measured 
3s3p( l Pi)nl series limit partly to the nl depend ence of the field - ioniza tion process in the 



experiment. Here we assumed a sharp n cutoff. ISchippers et al.l (120011 ) discuss the effects 
of a more correct treatment of the field-ionization process in TSR. Their formalism uses the 
hydrogenic approximation to take into account the radiative lifetime of the Rydberg level n 
into which the initially free electron is captured. 

Our theoretical calculations indicate there are no DR resonances due to 2 — > 3 or 3 — > 4 
core excitations below 44 eV, significant or insignificant. The two weak peaks above the 
3s3p( 1 Pi)n/ series limit at 43.63 eV are attributed to AN=1 resonances. 



6.2. Rate Coefficients 



The recommended rate coefficient of lArnaud fe Raymond! (119921 ) is in mixed agreement 
with our experimental results (Fig. [TU]) . For temperatures below 90 eV, their rate coefficient 
is in poor agreement. At temperatures where Fe XV is predicted to form in photoionzed gas, 
their data are a factor of 3 to o rders of magnitude smaller than our experimental results. At 
temperatures above 90 eV, the lArnaud fe Raymond! (119921 ) data are in good agreement with 
our combined experimental and theoretical rate coefficient. 
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As already implied by the work of iNetzer et al.l (120031 ) and iKraemer et al.l (120041 ) , the 
present result shows that the previously available theoretical DR rate coefficients for Fe XV 
are much too low at temperatures relevant for photoionized plasmas. Other storage ring 
measurements show similar differen ce with publ i shed recommended low temperature DR 
rate coefficients for Fe M-shell ions (lMullerlll999l ; ISchmidt et al.ll2006l ). The reason for this 
discrepancy is primarily because the earlier theoretical calculations were for high temperature 
plasmas and did not include the DR channels important for low temperatures plasmas. 

At temperatures relevant for the formation of Fe XV in photoionized gas, we find that 
the modified Fe XV rate coefficient of Netzer (2004) is up to an order of magnitude smaller 
than our experimental results. The modified rate coefficient of Kraemer et al. (2004) is a 
factor of over 3 times smaller. These rate coefficients were guesses meant to investigate the 
possibility that larger low temperature DR rate coefficients could explain the discrepancy 
between AGN observations and models. The initial results were suggestive that this is the 
case. Our work confirms that the previo usly recommended P R d ata are indeed t o o low but 
additionally shows that the estimates of Netzer et al. ( 20031) and Kraemer et al. ( 2004 ) are 
also still too low. A similar conclusion was reached by ISchmidt et al.l (120061 ) based on their 
measurement for Fe 13+ . Clearly new AG N modeling studies need to be carried out using our 
more accurate DR data (lBadneilll2006al ). 



Our state-of-the-art MCBP calculations are 37% lower than our experimental results at 
a temperature of 1 eV. This difference decreases roughly linearly with increasing temperature 
to ~ 25% at 2.5 eV. It is basically constant at ~ 23% up to 7 eV and then again nearly 
monotonically decreases to 19% at 15 eV. As discussed in §HJ a small part of these difference 
may be attributed to unsubtracted metastable 3 P contributions. But these contributions are 
< 10% at 2.5 eV, < 5% at 5 eV, < 2.0% at 10 eV, and < 1.4% above 15 eV (hence basically 
insignificant). Above 15 eV the difference decreases and at 23 eV and up the agreement is 
within < 10% with theory initially smaller than experiment but later greater. Part of the 
good agreement at these higher temperatures is due to our use of theory for the unmeasured 
DR contribution due to states with n > 80. 



7. Summary 

We have measured resonance strengths and energies for AiV=0 DR of Mg-like Fe XV 
forming Al-like Fe XIV for center-of-mass collision energies E cm from to 45 eV and compared 
our results with new MCBP calculations. We have generated an experimentally-derived 
plasma rate coefficient by convolving the measured MBRRC with a Maxwell-Boltzmann 
electron energy distribution. We have supplemented our measured MBRRC with MCBP cal- 
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culations to account for unmeasured DR into states which are field-ionized before detection. 
The resulting plas ma recombination rate coeffi cient has been compared to the recommended 
rate coefficient of lArnaud fe Raymond! (119921 ) and new calculations using a state-of-the-art 
MCBP theoretical method. We have considered the issues of metastable ions in our stored 
ion beam, enhanced recombination for collision energies near eV, and field-ionization of 
high Rydberg states in the storage ring bending magnets. 



As suggested by lNetzer et al.l (120031 ) and lKraemer et al.l (12004] ) . the present result shows 
that the previously available theoretical DR rate coefficients for Fe XV are much too low. 
Other storage ring measurements show similar difference s with published recommended low 
temperature DR rate coefficients for M-shell iron ions (jMullerlll999l ; ISchmidt et al.ll2006l ). 
We are now in the process of carrying out DR measurements for additional Fe M-shell ions. 
As these data become available we recommend that these experimentally-derived DR rate 
coefficients be incorporated into AGN spectral models in order to produce more reliable 
results. 

We gratefully acknowledge the excellent support by the MPI-K accelerator and TSR 
crews. CB, DVL, MS, and DWS were supported in part by the NASA Space Astrophysics 
Research Analysis program, the NASA Astronomy and Astrophysics Research and Analy- 
sis program, and the NASA Solar and Heliosperic Physics program. This work was also 
supported in part by the German research-funding agency DFG under contract no. Schi 
378/5. 
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Table 1. Energy levels for the n = 3 shell of Fe XV relative to the ground state. 



Level 


Rtiptetv (pV) 


3s3p( 3 P °) 


28.9927 


3s3p( 3 P 1 °) 


29.7141 


3s3p( 3 P 2 °) 


31.4697 


3s3p( 1 P 1 °) 


43.6314 


3p 2 ( 3 P ) 


68.7522 


3p\ l D 2 ) 


69.3816 


3p 2 ( 3 P) 


70.0017 


3p 2 ( 3 P 2 ) 


72.1344 


3p 2 (%) 


81.7833 


3s3d( 3 Pi) 


84.1570 


3s3d( 3 P> 2 ) 


84.2826 


3s3d( 3 P> 3 ) 


84.4848 


3s3d( 1 D 2 ) 


94.4875 


3p3d( 3 F°) 


115.087 


3p3rf( 3 P 3 °) 


116.313 


3p3rf( 3 P 4 °) 


117.743 


3p3d( l D° 2 ) 


117.601 


3p3d( 3 D°) 


121.860 


3p3d( 3 D°) 


123.346 


3p3d( 3 D°) 


123.565 


3p3d( 3 P 2 °) 


121.940 


3p3d( 3 P °) 


123.474 


3p3d( 3 P°) 


123.518 


3p3d( 1 F°) 


131.7351 


3p3d( 1 P?) 


133.2690 


3rf 2 ( 3 P 2 ) 


169.8994 


3rf 2 ( 3 P 3 ) 


170.1106 


3rf 2 ( 3 P 4 ) 


170.3612 


3d 2 ( 1 D 2 ) 


173.8992 


3d\ 1 G 4 ) 


174.4529 


3rf 2 ( 3 P ) 


174.2613 b 
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Table 1 — Continued 



Level 


Energy (eV) a 


3rf 2 ( 3 Px) 


174.3433 b 


3rf 2 ( 3 P 2 ) 


174.5416 


3d 2 ( 1 5 ) 


184.3712 



a Ralchenko et al. (2006) 
unless otherwise noted. 

b Churilov et al. (1989) 
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Table 2. Measured resonance energies Ei and strengths Si for Fe XV forming Fe XIV via 
N = 3 — > N' = 3 DR for .E cm < 0.95. Fitting errors are presented at a 90% confidence level. 



Peak Number E { (eV) S t (10" 21 cm 2 eV) 



1 (6.74 ± 0.05)E-3 189430.0 ± 20635.3 



2 


0.0098 


± 





,0008 


10078.0 ± 483.1 


3 


0.0196 


± 


0. 


.0008 


613.1 ± 56.8 


4 


0.0254 


± 





,0003 


743.9 ± 51.8 


5 


0.0444 


± 


0. 


,0002 


686.3 ± 37.9 


6 


0.0610 


± 





,0002 


2949.3 ± 39.0 


7 


0.1098 


± 


0. 


,0002 


805.5 ± 699.5 


8 


0.1674 


± 


0. 


.0014 


2424.3 ± 954.1 


9 


0.1943 


± 





.0018 


4408.5 ± 1213.1 


10 


0.2143 


± 


0. 


.0022 


4735.5 ± 750.9 


11 


0.2436 


± 





,0003 


4257.6 ± 132.6 


12 


0.2660 


± 


0. 


.0006 


4169.1 ± 339.0 


13 


0.2895 


± 





,0122 


213.9 ± 218.4 


14 


0.3102 


± 


0. 


,0074 


292.5 ± 188.6 


15 


0.3346 


± 





,0008 


1158.1 ± 118.6 


16 


0.3596 


± 


0. 


,0010 


943.5 ± 100.3 


17 


0.4154 


± 


0. 


.0149 


193.3 ± 230.2 


18 


0.4536 


± 





.0005 


8013.6 ± 328.0 


19 


0.4781 


± 


0. 


.0072 


706.9 ± 310.2 


20 


0.4988 


± 


0. 


.0072 


781.3 ± 303.5 


21 


0.5199 


± 


0. 


.0266 


216.7 ± 285.6 


22 


0.5433 


± 





.0290 


121.8 ± 270.4 


23 


0.6164 


± 


0. 


,0078 


136.2 ± 106.9 


24 


0.6599 


± 





,0006 


1269.1 ± 97.8 


25 


0.6992 


± 


0. 


,0010 


3090.3 ± 99.5 


26 


0.7385 


± 





,0010 


2068.5 ± 113.4 


27 


0.7943 


± 





.0006 


1594.4 ± 83.7 


28 


0.8406 


± 


0. 


.0006 


1740.6 ± 83.6 


29 


0.8830 


± 


0. 


.0006 


2164.2 ± 89.9 


30 


0.9232 


± 


0. 


.0013 


1420.7 ± 86.9 
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Table 3. Fit parameters for the total experimentally-derived DR rate coefficient for Fe XV 
forming Fe XIV via N = 3 — > N' = 3 core excitation channels and including the theoretical 
estimate for capture into n > 80 (n max = 1000). See § 4 for an explanation of the columns labeled 
"Experiment (I)" and "Experiment (II)". Also given are the fit parameters for our calculated 
MCBP results (n max = 1000). The units below are cm 3 s _1 K L5 for a and eV for 



Parameter Experiment (I) Experiment (II) MCBP 



Cl 


1 nvi7 /i 

l.U ltj-4: 


1 nvi7 A 
l.U 


i nvi7 /i 


C2 


Q Q(517 (5 

o.zoriv-o 


O Q(5T7 (5 

o.zoriv-o 


*7 1 Q17 O 


C3 


1 nni7 (5 
l.UUriv-0 


1 nni7 (5 


O (5*717 O 


C 4 


1 A (517 (5 

1.4ori/-0 


1 /I (517 £ 


Q 1 c;i7 o 


C5 


o tit? c 
Z. ( ihi-b 


o V7T7 (5 


1 (5017 1 


C6 


l.OlHv-0 


1 £ 1 17 (5 

l.olrlv-o 


r 0*717 /( 

O.o /rL/-4 


C7 


o nni 1 (5 
Z.yUriv-0 


o oni7 (5 
o.zyriv-o 




C8 


z.ooriv-o 


1 (5Q17 A 




C9 


k (5017 k 


A 1 AT? A 

4.14r!v-4 




ClO 


(5 (5(517 £ 

o.ooriv-o 


O 1 717 Q 
Z.Li Fj-O 




Cu 


(5 Q T 17 K 

o.olrlv-o 


(5 A H17 




Cl2 


V QQ17 K 


4.yoriv-z 




Cl3 


/i n*7T7 (5 
4.U / Hv-o 






C14 


c 0(517' (5 

o.yoriv-o 






Cl5 


z.o4ri/-o 






Cl6 


O OQ17 c; 
Z.Zocj-d 






Cl7 


5.27E-6 






Cl8 


2.40E-4 






Cl9 


2.22E-5 






C20 


2.56E-5 






C21 


7.40E-6 






C23 


4.35E-6 






C23 


5.51E-6 






C24 


5.50E-5 






C25 


1.42E-4 






C26 


1.00E-4 






C27 


8.32E-5 






C28 


9.61E-5 






C29 


1.25E-4 






C30 


8.61E-5 






C 3 1 


1.02E-4 
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Table 3 — Continued 



Parameter Experiment (I) Experiment (II) MCBP 



C32 


5.46E-1 


- 


- 


C33 


2.91E-3 


- 


- 


C34 


4.83E-3 


- 


- 


C35 


4.86E-2 


- 


- 


C36 


1.51E-1 


- 


- 


Ei 


6.74E-3 


6.74E-3 


4.12E-1 


E 2 


9.80E-3 


9.80E-3 


2.06E+0 


E 3 


1.97E-2 


1.97E-2 


1.03E+1 


E 4 


2.54E-2 


2.54E-2 


2.20E+1 


E 5 


4.45E-2 


4.45E-2 


4.22E+1 


E6 


6.10E-2 


6.10E-2 


3.41E+3 


E 7 


1.10E-1 


1.10E-1 


- 


E 8 


1.67E-1 


1.91E-1 


- 


Eg 


1.94E-1 


3.33E-1 


- 


E10 


2.14E-1 


9.63E-1 


- 


En 


2.44E-1 


2.47E+0 


- 


E12 


2.66E-1 


1.08E+1 


- 


E13 


2.90E-1 


3.83E+1 


- 


E14 


3.10E-1 






E15 


3.35E-1 






E16 


3.60E-1 






E 17 


4.15E-1 






Ei8 


4.54E-1 






E19 


4.78E-1 






E20 


4.99E-1 






E21 


5.20E-1 






E22 


5.43E-1 






E23 


6.16E-1 






E24 


6.60E-1 






E25 


6.99E-1 
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Table 3 — Continued 



Parameter Experiment (I) Experiment (II) MCBP 



E26 


7.39E-1 


E27 


7.94E-1 


E28 


8.41E-1 


E29 


8.83E-1 


E30 


9.23E-1 


E31 


1.00E+0 


E32 


1.16E+0 


E33 


1.62E+0 


E34 


3.14E+0 


E35 


1.08E+1 


E36 


3.82E+1 
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Fig. 1. — Fe XV to Fe XIV 3 — > 3 DR resonance structure versus center-of-mass energy E cra 
from to 1 eV. The solid curve represents the measured rate coefficient (av) which is the summed 
DR plus radiative recombination (RR) cross sections times the relative velocity convolved with the 
experimental energy spread, i.e., a merged beam recombination rate coefficient (MBRRC). The 
dotted curve shows our calculated multiconfiguration Breit-Pauli (MCBP) results (n max = 80) for 
ground state Fe XV (top plot) and 3 Pq metastable state Fe XV multiplied by a factor of 0.06 
to account for the estimated 6% population in our ion beam (bottom plot). To these results we 
have added the convolved, non-resonant RR contribution obtained from semi-classical calculations 
Jschiooers et al.l B. The inset shows our results for E cm from 5 x 10- to 1 x 10- eV. 
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Fig. 7. — Same as Fig. 2 but for E cm from 23 to 36 eV. The dotted curve shows our calculated 
MCBP results and the thin solid curve shows our calculated MCBP results reduced by a 
factor of 1.31. 
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Fig. 8. — Same as Fig. 7 but for E cm from 35 to 45 eV. The weak resonances above 44 eV 
are attributed to AN=1 DR. These are not included in either our experimentally-derived or 
theoretical Maxwellian rate coefficients. 



33 



I I I I I — I — I— 



I I I I I — I — I- 




> 

3. 

<? (1) 

LJJ 

w 
w 

CO 



< 



CD 

-I— ' 

O CD 

^ o 



III — I — I 1 1 — 



I— I— I — I — I — I — I 1 1— 



I I I I I — I — I — I — 



( t s e wo 8 _(H)< AD > 



Fig. 9. — Measured and fitted Fe XV to Fe XIV 3^3 resonance structure below 0.07 eV. 
The experimental MBRRC results are shown by the filled circles. The vertical error bars 
show the statistical uncertainty of the data points. The solid curve is the fit to the data 
using our calculated RR rate coefficient (dashed curve) and taking into account all resolved 
DR resonances. The dotted curves show the fitted DR resonances. At E cm = 0.005 meV the 
difference between the model spectrum a and the data is 1 + (A«/a ) = 2.5. 
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Fig. 10. — Maxwellian- averaged 3^3 DR rate coefficients for Fe XV forming Fe XIV. The 
solid curve represent our experiment ally- derived rate coefficient plus the theoretical estimate 
for unmeasured contributions due to capture into states with n > 80. The error bars show 
our estimated total experimental uncertainty of ±29% (at a 90% confidence level). No error 
bars are shown below 1 eV for reasons discussed in § HI The thin solid curve represents our 
experimentally-derived rate coefficient without the two lowest energy resonances included. 
The dash-dotted curve represents our experimentally-derived rate coefficient alone (n max = 
80). Also shown is the recommended DR rate coefficient of Arnaud & Raymond (1992; thick 
dash- dot- dotted curve) and its modification by Netzer (2004; thin dash- do t- dotted curve). 



The fi lled pentagon at 5.2 eV represents the estimated rate coefficient from iKraemer et al. 



( 12004] ) . The dashed curve shows our MCBP calculations for n max = 1000. As a reference 
we show the recommended RR rate coefficient of Arnaud & Raymond (1992; dotted curve). 
Neither the experimental nor theoretical DR rate coefficients include RR. The horizontal 
l ine shows the temperature range over which Fe XV is predicted to form in photoionized gas 
f kallman fc BautistalEoOlh . 
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